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Abstract 

Using dispersion theoretic techniques, we consider coherent long range forces 
arising from double pseudoscalar exchange among fermions. We find that 
Yukawa type coupling leads to spin independent attractive potentials 
whereas derivative coupling renders spin independent repulsive poten- 
tials. 
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I. INTRODUCTION 



Many extensions of the Standard model predict the existence of hght scalar particles. 
The axion may be the most debated one but there are also approximate Nambu-Goldstone 
fields associated with family symmetries, or moduli fields, or dilatons, or superpartners of 
the gravitino Exchange of such particles by ordinary matter will induce forces that 

extend over the Compton wavelength of the particle 0,^,0. However, the effect will be 
felt by bulk matter only if the potential is spin-independent so that forces can add up co- 
herently over macroscopic distances. Now, a pseudoscalar particle, such as the axion, is 
coupled to fermions via a 75 which, in the nonrelativistic limit, flips the spin. Therefore, 
single pseudoscalar exchange leads to spin-dependent forces that do not extend over macro- 
scopic unpolarized bodies A double exchange of pseudoscalars on the other hand can 
coherently sum over a macroscopic sample of matter because it can leave the spin unfiipped. 
The explicit form of these forces has been derived and their phenomenological consequences 
explored in previous work in the context of nonrelativistic "old fashioned perturbation the- 
ory" 10]. Here we reopen the question of pseudoscalar mediated forces in the light of the 
powerful dispersion theoretical techniques devised by Feinberg and Sucher and collaborators 
that make extensive use of full relativistic quantum field theory 1^,0. 

In section ^ we give the necessary theoretical background which is nothing but a brief 
resume of the seminal work by Feinberg and Sucher. Section |TTI| is devoted to the Yukawa 
type interaction and section deals with derivatively coupled scalars. We shall see that 
the two interactions produce quite different potentials |jTO|, a fact that could not have been 
derived in a purely non-relativistic framework |[Tl|. The paper ends with a brief summary 
and conclusions contained in section 0. 
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II. DISPERSIVE FORCES 



Following the general strategy devised by Feinberg and Sucher , we define a potential 
in a given Quantum Field Theory by equating the scattering amplitude for a two body 
process that follows from the usual Feynman rules, with the transition amplitude associated 
to a Schrodinger type equation solved a la Lippmann-Schwinger. Let us be exphcit and 
consider elastic scattering of particles A and B with four momenta Pa and in the initial 
state and p'^ and p'^ in the final state. The Mandelstam variables are then: 

S = {Pa+ Pbf t = Q'^ u = {pa- p'bf (1) 

with Q = p^-p'^ = -pt,+ pI 

In the CM. we write the momenta as 

Pa = {Ea, p) Pb = {Eb, -p) 
p'a = iEa,p') p', = iEb,-p'). 
Now s = W"^ where W = Ea + E^ and t = — with Q = (0, Q), whose physical region 

is 

s > So and -Ap'^ <t<0 (3) 

where 

So = {rua + rribf p^ = [s - {rria + rribf] [s - {nia - m^f] j 4s. (4) 

The transition from initial state i to final state / is described in Quantum Field Theory 
by the transition matrix element 

Ts. = NjMf.N, (5) 

where Nf^i are normalization factors of one particle states and Ai fi is the invariant Feynman 
amplitude. 

The definition of our potential follows now from identifying this transition amplitude 
with 
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T/. = (p', -p'\V + V{W-ho-V + ze)-' V |p, -p) ^ (6) 

where ho is the sum of the free Dirac Hamiltonians for particles A and B. Here all quantities 
are referred to the CM. 

The Feynman amplitude Ai is understood as a series expansion in (even) powers of the 
coupling constant associated to single, double, . . . particle exchange in the t-channel. We 
assume that the potential V also admits a series expansion 

So we determine V order by order in perturbation theory through 

(p',-p'|02)|p,-p) =M^ (8) 



(P, -P'l V^'^ |p, -p) = Mf^ - (p', -p'l F(^) {W-ho + le)-' V^^^ |p, -p) . (9) 



,(4) 

Here we should point out a technicality. The potential we are after must have the 
form [jr^l 

l^ = A++[/A++ (10) 

where A+_|_ = A+.^A+.b is an operator that projects on the positive energy states of h^. 
Because A++ | p, — p) =| p, — p) and A^_^ = A++ we can rewrite the previous equations that 
determine the potential as, 

(p',-p'|f/(2)|p,-p) =A^(? (11) 



(P, -p'l C'*" |P, -P> = J^t! - (p', -p'l C'" (W - li« + A++ C/<'' |p, -P) (12) 



In principle the above equations permit an iterative determination of the potential to 
the desired order. However we are not done yet because we would like to have our potential 
in position space and what we have is the operator U in the momentum representation. 
Therefore, we should Fourier transform our results back to configuration space, i.e. we wish 
to find U^^Hr) such that 



(p', -p'l |p, -p) = j dr e'^ ' U^''\r). (13) 

Inversion of the above equation, however, requires knowing the function Ai for all values 
of three momentum. But we only know the scattering amplitude on shell, i.e. for p^ = p'^. 
We can use the fact that A4{s,t) is an analytic function of t and so analytically extend its 
domain beyond the physical region, i.e. for all values of = t. 

Suppose A4{s,t) is analytic everywhere except for branch cuts on the real axis and 
furthermore, it vanishes for large \t\. Then, using Cauchy's theorem, we can write 

M^-^^'-Tdt'^^ (14) 



^Jto t'-t 

A.(^) ^ i r (15) 

where p{s,t) = -^^^ is the spectral density and [A4]^ is the discontinuity of j\4 across the 
cut. Only the piece of the amplitude arising from the right hand cut will be of interest to 
us for only this piece leads to a long range potential [|l^ . Assuming that the basic relations 



eqs ([TT|) and ([l2|) hold also in the extended domain, we can Fourier invert them as follow^, 

[/W(r; s) = ^3 [dQ e-Q-"- (s, -Q^). (16) 

(27r) J 

We use now the spectral representation given before to obtain 

1 f°° 

U^^\r;s) = —- / dt p^''\s,t) e-^'"- (17) 

where, to reach this final form, we conveniently changed the order in which integrals were 
done. In short, obtaining long range potentials amounts to calculating t-channel disconti- 
nuities in Feynman diagrams and performing a Laplace transform. We shall see how things 
work out in detail as we do our specific calculations in the next two sections. 



^Notice that our generalized potential will depend on the parameter s. 
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III. THE YUKAWA COUPLING 



Our starting point is the Lagrangian density, 

£L = -z^7^(x)7'v^(a;)$(x) (18) 

where ■?/' is a fermion field and is the pseudoscalar field which we take to be masslessQ. 
The potential associated to single particle exchange is easily obtained from the discontinuity 
associated to the diagram in Fig. 1. The spectral density function p^'^\s,t) is in this case 

p^'\s,t) = 7rg'u{p'J^'uipa) uipi)j'uip,) S{t). (19) 

After Laplace transforming we get the relativistic potential operator 

U^'' = ihl ihl (20) 

where subindices make explicit that Dirac matrices act either on spinor A or spinor B. 

The nonrelativistic limit of equation above leads to the well known spin-dependent po- 
tential 

= 4. r(2,t)(2,n.) "'^ ® "^^ 

with ma and m;,, the masses of particles A and B respectively. What we are really inter- 
ested in is U^^\ i.e. the potential due to two-particle exchange. To this end we need the 
discontinuities of diagrams in Fig. 2 and the discontinuity of the subtraction term in eq (p!2D 
(iteration of the lowest order potential t/*-^-*). 

The Feynman amplitude associated to Fig. 2 can be written 



2!(27r)4 J ' k^ + iek'^ H 



M'::{-k,k';Pa)M^'{k,-k';P,) (22) 



^For scalars with mass, the long distance potentials are damped with Yukawa exponentials. In 
this case, our results are valid for distances on the order or smaller than the Compton wavelength 
of the exchanged particles. 



in terms of the Compton amplitude, depicted in Fig. 3 

[k, k'- P) = / u{p' 



+ 



u{p). 



(23) 



2p-k 2p-k'_ 

Making use of the Dirac equation and trading the propagators for their discontinuities, 



I.e. 



1 



k"^ + it 



-2m 5 (k^) e (A;°) 



(24) 



we arrive at 



87r2 



Pa ■ jk' - k) ^ 
2pa-k' Pa-k 



■ jk' - k) ^ 

2pb- k' pb- k 



Ub 



(25) 



with the two particle phase space explicitly given by 



d^ = 6{Q-k-k')6 {e) 6 (A;'') 6 (A;°) 6 (A;'°) d^k d^k'. 



(26) 



It is convenient to do the integrals in the CM. of the pseudoscalars, i.e. to go to the 
t-channel, and then use crossing symmetry to recover the original amplitude. We follow here 
the notation in reference |T^ where they deal with a related problem. First define momenta 



as, 



Pa = 

Pb = -Pb 
k = 



Vt 1^ 

2 



Pa = -Pa 

P'b = 
k' = 



Vt _i 

2 ' * 



(27) 



introduce next the unit imaginary vectors 



p = i iarria p 
p ' = 2 ^fcmf, p' 



(28) 
(29) 



with ^ab = — A ^iid p, p' are unitary complex vectors verifying p ■ p = — 1 so that 
all particles are on-shell. 

Now the discontinuity can be put in the form 
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where, to simphfy expressions, we use 

b = rua^a rrib^b 
Xa = P ■ k Xb = p' ■ k 

= 2«„,tl„.i, ■ (31) 

The integration to be carried out is an angular average. We use the shorthand: J ^/ = 
(/). Hence, the discontinuity takes the form 

with ri"'=(^^^ k'^k-'Y 

Lorentz covariance dictates the following decomposition 

T^^'' = a, PjiP: + a^P^PI^ + a, {P^Pi; + P.^P^) + a,g^^ 

+ a^Q^Q'' + a, {Q^P^ + P^Q^ + a, {Q^P^ + P^Q^) (33) 

in terms of the three independent momenta 

Pa=Pa+p'a Pb = Pb + p[ Q = k + k' . (34) 

The coefficients Oj can be found to be combinations of scalar integrals as shown in the 
appendix. Now, in the CM. of the incident particles, the relations 

u'aPbUa u'bPbUb ^ 4: Ulb < {W^fo - Ula) Ua u'^Ub 
u'aPaUa u'bfbUb = 4 m„mfe u'^Ua u[Ub 

u'aPbUa u'^PaUb = 4 u'^ {W^fo - ma) Ua u'^ (H^7o - rUb) Ub 
u'aQUa = 

u'bQub = (35) 



are easily established with the help of the Dirac equation. This leads directly to 
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1 CI 

[M'^^^]^ = -^-g^ u'^u[ [4m„mb (2a3 - - aa) + Ami^W^^ (as - ag) 

+ 4m, 1^7° (ai - a^) + 7°7° {AW^a^^ + 04) - 7a7fe «4] (36) 

This discontinuity is a complex function since the are complex and hence adds an 
imaginary component to the spectral density which would finally contribute an imaginary 
piece to the potential. Inspection of equation (^) immediately tells us that the offending 
piece comes from the imaginary parts of the Oj. But we should recall that we still have 
to subtract the contribution from the iterated lowest order potential. It turns out that 
its imaginary part exactly cancels the unwanted contribution coming from equation (p6|). 
Indeed we have explicitly checked this to be the case. However, in order to make this paper 
not too lengthy, we do not include the intermediate steps of the calculation. We only report 
on the result, i.e. 



p2 



(4p2 + tf 



{Eat - ^a) {E,t - m,) 



UaUb- (37) 



The relevant contribution to the long range potential comes from the real parts of the a,, 
that is the imaginary part of (^), once the contribution of the iterated potential has been 
subtracted. Let us elaborate on the iteration amplitude. 

Mi = (p', -p'l f/(2) iW-ho + te)-' A++ U^^^ |p, -p) . (38) 

This formal expression can be recast in the explicit forniQ 

^/ = 1^ / I'dl I ^ ul'ul' 7° (^7° -7a-l-ma) 

7° + 76-1-^6) UaUk Cip, l)-^^ (39) 



^This integral as it stands is infrared divergent. A fictitious mass regulator is understood to be 
introduced in the scalar propagators which is set to zero after the integrations are performed. 
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with q = p - 1 and q' = p' - 1 and C{p, I) = j;^j;/(vK-vK'+ie) • 

The integration over momentum 1 reflects the fact that we have inserted a complete set 
of plane wave intermediate states in (^). We have also used 

A«(l) = ^^ + ^°^^-^,' + -'>° (40) 

where E'^ = a/?7i|~+7^. 

The iteration amplitude can be conveniently put as follows 

4 r 

= 1^ y ^'dl y'A { {Kla - ma) {Ei^l - m,) /: - 7. • V {E[^l - m,) 

+ (e;7° - m,) 7, ■ V - 7h^T,, } Uau, C{p, I) (41) 



where 



' 47r q'2q2 ^ ^ 



V - / - 1 -r,-. (43) 



dVL ^ I I 
Atc q'^q" 

and they can be found in the appendix. 

A little bit of Dirac algebra and the results in the appendix allow us to write the discon- 
tinuity of TVI/ as 

,4 



[Mj], J fdl v!X { {Kll - m^) {Kit - m,) 



[(E,7a - ma) - m,) + (E^o - nia) {E,^', - m,)] 

A {Eall-ma) {E,^l-m,) 



Ap'^ + t 
4 fJp' + l'f 



The explicit form for [C]t is given in the appendix. Note that in (^Sf ) we have 

^ p( T,, ) -i7r^(Vr-iyO. (46) 



W-W' + te " \ W-W 
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The Dirac delta piece gives a contribution that, as aheady advertised, will exactly cancel 
the real part of the fourth order discontinuity function (pT]). The principal part integral can 
be cast in the form 



[MA, 



p 



dx 



x" 



C{p,l) U^u[AfUaUb 



(47) 



where we changed the integration variable via the relation 



a' + b'x 1 



(2p2 + t) + ^/t{Ap'^ + t)x 



(48) 



and we used the shorthand 



Af ^ (4p^ + tr [K^l - m,) (Eh,° - m,) - 2(4p2 + t)(/ + 

[(E,7° - rua) - m,) + (^^7° - m,) (^,7° - m,)] 

+4 (2(p2 + e f - l\Ap' + t)) [Ball - ma) {E,^l - m,) 
-la ■ 7.(V + t) ((4p2 + t)l' - {j? + l^f) . 

It is convenient now to split the function C(p, /) as 



c(p,o = Ci(p,/) + C2(p,/) 

2 



(49) 



1 (p^ + l'^ + ml + ml W 

7— If 7— 1/ T T 7- I T— I 7—1 7— 1 » 7— 1/ 



(50) 



l^(p2-/2) E'^E'.W \ EaEb + E',E', W + W'J' 

The integral above cannot be done exactly and we will expand the integral in a power 
series in t and p^. This is a licit procedure because we will perform a Laplace transform 
that heavily weighs the small t region of the spectral function when determining the long 
range (large r) potential and, we will eventually take the non relativistic limit of the potential, 
i.e. for p^ ~ 0. Furthermore, each extra power of t or p^ implies a correction to the potential 
with an extra power of r~^. We see from equation (|50|) that the calculation will involve 
doing integrals of the type 



Iclijl) 



dx 



I'' 



1 Vl^P" - 



(51) 



11 



Ic2(n) 



dx 



r. 



(52) 



-1 Vl - X 

The explicit results of the integrals needed in our calculation are also given in the ap- 
pendix. Armed with all this artillery we find for the Ci piece of the discontinuity [A^/]*: 



8Vi{Ap^ + ty{ma + mb) 



mairib 



4t- 



pH 



1 



1 



+ 7° maUib 



+ 



mpub 

t^p^ —ml — mlnib + maml + 2m^ 4 
~t~ ~ 7 ~t~ ~ ~^ q tp 



32mi Ami 



2m^ml 



+ 



_^^^ 2K-m.) ^^, 
rriamf 

, 2m^ + mlmb - mgmf - mf ^ 4 



32ml Ami 



2mlml 



+ lllb ^atrib 



4^ + — + — 
1 1 



32m^ 32m^ Smlmf 
2^2 



J.2 2 



1 1 1 



4m^ 4ml 



m'^ml 



-tp 



4 2m^ + m^mb + mimi + mam^ + 2m^ 



and for the C2 piece: 



2tp' + -- 



2m^ml 
tp^ 



2 Amamb mamb 



UaUb 



(53) 



7a ruamb 



16"\/t(4p2 + ty{ma + mb)mamb 
ml + mamb + ml^^^^^ 



-I -I r 



IQmlmb 
ml + mamb + ^^^^^2 ^ ^^2 



+ 7^ m„m6 

ml + mlmb + 2m^m^ + mamf + 



7°7.° 



la' lb 



IQm^ml 
ml + m„m6 + ^"|^^4^2 ^ ^^2' 



i(4p2 + i)2 



(54) 



16mamb 

In both equations above we kept only up to the powers of t and that will be needed, 
either in this section or in the next section, to obtain the leading two-particle exchange 
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potential. In this respect equation (|54D does not contribute to the leading potential just 
under scrutiny. To obtain the spectral density p*^'*-' , we must finally perform the subtraction 



(55) 



Recall that what enters p^^^ is the real part of the a, in (PB|). Although the integrals 
that go in the Oj are exactly given in the appendix, the required subtraction ( ^5|) and final 
Laplace transformation ([TtD demand that we here also expand the integrands in a power 
series in t and p^. After some lengthy algebra, we arrive at the final form for the imaginary 
part of [M]t 



Q [M], = Q [M]f'' + $5 [M] 



(56) 



where 



53 



odd 



+ 7° 



4(m„ + mb)yt(4p2 + t)^ 



2t 



-3m^ + 2m^m^ 



rriamb 



^tp^ 



— 3m^ — m\mb + 



mlml 



rriaml 



3m^ 



+ 



-3m^ — m^nib + 



•im^rril 



rriaml 



-2t + 



rria 



mlrrib 



nib 2 , + 4mb 4 
tp +— —tp 



3m^ + bm\mb + ^a^^ 



Arriamf 
2 mgml - m^ ^2 2 



+ 



3m^ + 5mlmb + m'lm'l — mam\ — 3 



GArriaml 



-2t + 



Arria + 3mb 



mb-rria. 2 , 

o-^P ^ — 2 

mam^ Arriamb 



tp' 



—m\ — rriamb + Tn^ml + bniaml + Smf ^ 2 



+ 7°7.° 



2t + 



— rriamb + rriaml + dmamf + 3mf^^ 

tp4 
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o 4 4 ^ ^ 



+ la-lb 



+ 



, + AmaTrib + ml ^ 
H ^ , 



marub 4:mamb 
3ml + bmanib + 3ml ^2^2 , 2m^ + mamf, + 2m? 3 



+ 



32m3 



(57) 



contains the odd powers of and 

^ i =^«a%< a 2 + 1 

^ 47r bmpub bmam^ Aniamb 



mqiTib + {{mg + mb)7° - m^) {{mg + mb)7g - m^) ] 



12m^m^ 



(58) 



contains the even powers of y/t. We did this separation to emphasize that, after the subtrac- 
tion in equation (^), only the term (^) survives to leading non- vanishing order. Indeed, 
equation (|57|) coincides exactly with equation (|53|) if we neglect terms beyond or tp^. 

The final step involves the Laplace transformation indicated by equation ([T7D . Using the 
general formula 



»,2rt+2 



(59) 



we get 



9" ^1 ^1 J {mg + mb)-il - "^b ("^a + "^6)7° - , 7a7|^ 



12m^m^ 



Gmginl 



Arrigmb 

UgUb (60) 



which leads, in the non relativistic limit and concentrating only on the spin-independent 
terms of (0), to the long-range attractive potential: 

,4 
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QAn^r^niamb 

where this operator is supposed to act between two-component Pauli spinors. 



(61) 
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IV. THE DERIVATIVE COUPLING 



In this section we consider the interaction Lagrangian 



(62) 



which is how Goldstone bosons couple to fermions. 

This derivative couphng leads to the same one particle exchange Feynman amplitude as 
before and therefore to the same lowest order potential ( pOD . Hence the iteration amplitude 
will be also identical. However, the two particle exchange amplitude (see, Fig. 2) is different 



because the Compton amplitude that goes into ( p2D is different. Indeed, the Compton 
amplitude, corresponding to Fig. 3, is now. 



{k, k'- P) = / u{p') 



+ 



u{p). 



(63) 



2p ■ k 2p ■ k' m 
This amplitude differs from (^) by an extra term proportional to m~^. 
We introduce this amplitude in (^) and replace the massless propagators by Dirac deltas 
to obtain the discontinuity function: 



87r2 



u' 



Pa-ik'-k)^ 
'^Pa -k' Pa - k 
1 



1 



Ub 



Pb ■ jk' - k)f 
2pb- k' Pb - k_ 

Pb - {k' - k) 



niamb rua 



2pb- kpb- k' 



Ub 



1 



Vlb 



Pa ■ jk' - k) ^ 
2pa-kpa-k'' 



UaUf^Ub 



(64) 



The first piece is exactly what we had in the last section. We call [AAl'^'^^j^ the extra 
added piece that involves the integrations 

1 



Stt^ 



rriamb 



^9 

IGnniamb 



U^Ua U^Ub, 



(65) 



and 



a 



d^ u'aUa U^jtUb 



Pb - jk' - k) 
2pb- kph- k' 

iVt^brribXb 



I (I + i^bmbXbr) 



Ub + (a ^ b) 

+ (a^b) 
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A . /I 



1 — Tfe arctan I — '^I'^a ^'^,^6 + (a 6) . (66) 



The last line in equation ( pB]) is reached by demanding Lorentz covariance to write 

{2 ^ ^ ^ ^^^^ 
by solving for the coefficients as explained in the appendix, and by using the Dirac equation. 
Putting things together, 

[AA<(^)] = u'^Ua u'b [4 - arctan ( -) ) 

1 / f l\\ 1 , , 

UaUb (68) 



+ ^ 1^1 - r;, arctan ^— ) ) - 1 



To leading order in t and we have 



In the nonrelativistic limit this contributes the quantity 



,4 



A\/i,^) = -i- / Ape-v^^rft=— -4 (70) 



which exactly cancels contribution (|6T1), i.e. 

C.. = + ^Kl'^ = + ^ • (71) 

Hence, to find the form for the potential in the case under scrutiny, we must go to the 
next order in our series expansions. What we need now is to consistently take into account 
the previously neglected terms in the spectral density 



(4) _ ^ 

^der - 22 



(72) 



So we collect the relevant pieces in (|57D , (^31), and (0) in addition to 

1(4) ^5''^V^ ^1:^1 1 , 1 
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which is the next to leading term in the expansion of (|68|). The result is 



(4) 
Pder 



M^f^ + /\M - Mi 



2i 

g^mambVt 

—, r u„Uu \ 



+7; 



— 3m^ — 2m\mh + brriaml + 3m^ 



3"^a + '^'n^a^b + 7m^m^ + 2mami + 3m5 
+spin-dependent terms of the kind 7^ • 7^} UaWb- 



(74) 



This spectral density nonetheless gives a vanishing spin-independent potential in the 
static approximation, i.e. 



(4) 



der:nr 



+ (r-^) . 



(75) 



The first non-vanishing contribution to the spin-independent potential arises from the 
part in the spectral density which is linear in t. Indeed, the explicit form of the spectral 
density reads: 



P 



(4) 



2i 



3_ 1 
Stt 240m4m^ 



4 [(— 3m^ — 2rn?jn}j + 2m^m^ — 2mam\ — 3m^) 
6m^ + Sm^mb + Am^mb + 3mam^ + 6m^j tj 
+ 7° 2 (ma + mb) [{6ml - 2mlmb + Sniaml - Ami) 

+ [—12ml + Sm^mf, — Amaml + 3ml) t] 
-f 7° 2(ma + mf,) [(— 4mj^ + Sm^mb — 2maml + 6ml) 
+ i^ml — Am^mb + 6maml — 12ml) ^] 
+ 7a7° [(8"^a + 2m||mfe - 12m^m^ + 2mamj| + 8m^) 

+ (-6m^ + 2mlmb - 9m^m^ + 2maml - 6ml) t\ } UaWf, 
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(76) 



where no iterated second order amplitude contributes to this order, and where we picked 
the term proportional to t in the expansion of ( |65D and (p^). 
If we use now 

< 7° ^ < (1 + Oip')) Ua (77) 

and pass to the static limit, we find 

"-=324?'"®'' 
which, upon Laplace transformation, leads to 

for the desired spin-independent long range potential. Note that, as opposed to the Yukawa 
type coupling potential (pT|), the derivative interaction leads to a repulsive potential. 



V. CONCLUSIONS 

Very light particles can mediate forces extending over distances on the order of their 
Compton wavelength. If this range is macroscopic, unpolarized bulk matter will only expe- 
rience the effect of spin-independent interactions. It is a well known fact that the Yukawa 
potential due to pseudoscalar exchange depends on spin and as a consequence no coherent 
effects do arise on a macroscopic scale, unless of course our sample is polarized ||5|. How- 
ever, residual (Van der Waals type) forces may arise between macroscopic bodies in the 
case of pseudoscalar mediated interactions, due to the exchange of two quanta " at the same 
time". The double helicity flip involved eventually makes the resulting effective potential 
spin-independent |]^. 

In the preceding sections we have established, with the help of the formalism developed by 
Feinberg and Sucher, the large distance behavior of such residual forces, i.e. those associated 
to double pseudoscalar exchange. We have considered two different basic couplings of those 
scalars to matter fermions. On the one hand we took the ordinary Yukawa coupling (e.g., 
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this is the way the Higgs particle couples to fermions) and on the other we considered the 
derivative coupling (e.g., the axion-fermion interaction). Both couplings reduce to the same 
spin-flip interaction in the static non-relativistic limit. And both interactions also produce 
identical one-particle exchange (spin- dependent) potentials. In spite of this fact, we have 
explicitly shown that the spin-independent two-particle exchange potential is substantially 
different in both cases. Indeed, for Yukawa coupling we derive a attractive long distance 
behavior whereas for the derivative coupling the potential, that is now repulsive, falls off as 
Since, as emphasized, double exchange will lead to residual macroscopic effects, these 
effects will be quite different in both cases. So, we have found still another instance where 
the interaction of pseudoscalars to fermions can be discriminated. Other places are, soft pion 
emission in proton-proton scattering or axion bremsstrahlung in a supernova core P,p!5[|. 

Of course, the effects just reported are extremely small for the light scalars presently 
contemplated in particle physics, such as the axion, and thus their experimental detection 
is beyond reach of present technology. However, there is much activity and interest on 
the experimental front and experiments are designed and performed that explore the sub- 
centimeter and sub-millimeter regime with an ever increasing sensitivity P^JT7[ ]. And, on 
the theoretical side, the completion of the Particle Physics Paradigm may still require new 
superlight scalar particles to exist. 
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APPENDIX: 
The coefficients 

The coefficients in the tensor decomposition (^) read: 
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^ t ih - 2/5 -h + Ayh - y^h - y^h) 

i6miea{i-y'y 

_ t{h-h- 2/6 + Ayh - y'h - y^h) 

lQmie,{l-y'? 

^ t {-h - yh + 2yh + 2yh - ^y^h + Z/'/i) 
166(1 -y2)2 

t{I,-I,-I, + 2yh-y%) 
«^ = 4(W^ 

with y = p ■ p', in terms of various angular integrals in the set 



dadb I \ dadb I \ dadh 

/ = / b \ J _ / -^g-^fe \ J _ / ■^g-^o 



/6=(5^). (A2) 

These results are obtained after repeated contraction of (^) with the independent mo- 
menta in our problem and after solving the resulting algebraic system of equations. 

The angular integrals Jj are given next. The explicit calculation of Jj for i < 5 is given 



m 



13 1 Is can be found in [14| and Jg follows trivially from Is by interchanging a and h. For 



the first two integrals from the set ([A2|) one gets 



where we have defined, 



/„ = ^±I^ /, = ^±I^ (A3) 



F+ = ±N'^ arctan f — | - A^~^ arctan f — 
i^Jp^s + byt 



b 

D± = y± Tan (A4) 



and the rest is given in terms of Jq and Ji by 



1 n 
— arctan — I — t^Iq 
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1 / 1 . 

i3 = — arctan I — ) — 



li = 1 — Ta arctan ( — ] — Tb arctan ( — ) + r^r^ Jq 



n 



-^5 = y ( 1 ~ ^6 arctan 



lQ = y[l-Ta arctan ( — ) ) - rlh. 



(A5) 



Angular integrals entering the iterated amplitude 



We start with (^21). Its discontinuity [L\t has been obtained in |T3[: 

1 



IT 



lVt^{ll-P){l^-ti) 



e {il - P) e (/2 - ii) 



with 



ll = a' = + t, b' = ^t{4p^ + t). 



(A6) 



(A7) 



For p = I , simphfies to 



(A8) 



We turn now to (^). It has been calculated in |T^. It is symmetric in p and p'. Hence 
it can be cast in the form 



V = t;P with P = p + p ' 



(A9) 



where v is given by 



(p + 0' 



(AlO) 



The associated discontinuity is then 



Finally let us discuss (0). This integral admits the general decomposition 



(All) 
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Tij = aP^Pj + hQiQj + c5ij (A12) 

in terms of the vectors P, defined above, and Q = p — p'. Clearly, we can write 

PiP^T^ = a{P^Y + cP^ 
Q.QjT' = b{Q^Y + cQ^ 
6ijT'^ = aP^ + bQ^ + 3c 



(A13) 



(we used: P ■ Q = 0). The scalar integrals on the left are, respectively 



PiPjT^ 



4^(q'2) (q2) 



{p - 0' 



4p2 



+ 



tin 



p. 



{P + 0' 

{p - ly 



(AM) 



Now the three equations can be solved for a, b, and c. Recall that we only need the 
discontinuity of (H), i.e. the discontinuities of a, b, and c. They are. 



+ t Ap'^ + t 



jp' + Py 
Ap'^ + t 



(A15) 



Integrals Id and Ic2 



Here we display the explicit solutions of ([51| ) and (|52D for n = 0, 2, 4, 6: 
/ci(0) = 

/cl(2) = -TT 



/el(4) = --(4p2+t) 

/ci(6) = -^((4p2 + t)^ + 2(2/ + tf) 



(A16) 
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Ic2{0) = TT 

7^2(6) = — (16p^ + 48pH + 30pH^ + 5t^). 
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Figure Captions 

FIG. 1 Lowest order scattering amplitude. Single pseudoescalar exchange. 
FIG. 2 Diagrams contributing to the 0{g'^) terms of the potential. 
FIG. 3 Compton scattering amplitude diagrams. 
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FIG. 1. 
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FIG. 3. 
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